Probing thermoelectric transport with cold atoms 
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We propose experimental protocols to reveal thermoelectric and thermal effects in the transport 
properties of ultracold fermionic atoms, using the two-terminal setup recently realized at ETH. We 
show in particular that, for two reservoirs having equal particle numbers but different temperatures 
initially, the observation of a transient particle number imbalance during equilibration is a direct 
evidence of thermoelectric (off-diagonal) transport coefficients. This is a time-dependent analogue 
of the Seebeck effect, and a corresponding analogue of the Peltier effect can be proposed. We reveal 
that in addition to the thermoelectric coupling of the constriction a thermoelectric coupling also 
arises due to the finite dilatation coefficient of the reservoirs. We present a theoretical analysis of 
the protocols, and assess their feasibility by estimating the corresponding temperature and particle 
number imbalances in realistic current experimental conditions. 
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Thermoelectricity has been a recurrent theme in con- 
densed matter physics. The increasing demand for sus- 
tainable energy sources, as well as progress in materials 
science, have triggered a marked increase of interest and 
research in the subject over the past two decades [1-4]. 
A better understanding of fundamental processes con- 
trolling thermal and thermoelectric transport is likely to 
bring progress in the field at large. 

Ultra-cold atomic gases offer a remarkably clean and 
controllable set-up to investigate interacting quantum 
systems. Phenomena involving the transport of atoms 
have been the focus of several experiments in this context, 
for example probing the reaction to external forces [5, 6] 
as observing Bloch oscillations in optical lattices [7, 8], 
or the expansion dynamics of atoms in disordered [9- 
11] and lattice potentials [12]. The transport of impuri- 
ties [13-15] and spin diffusion [16] have also been inves- 
tigated. Recently, increasing interest has been devoted 
to closer analogues of mesoscopic transport devices, e.g. 
modelling 'quantum pumps' and 'batteries' [17-19]. In a 
recent experiment transport of fermionic atoms between 
two reservoirs connected through a tunable constriction 
was realized [20]. 

In this letter, we go beyond atom transport and de- 
sign a proposal to observe both thermal and thermo- 
electric transport in the context of cold atoms, in the 
two terminal geometry of Ref. [20] (Fig. 1). We sug- 
gest protocols specifically geared at revealing the off- 
diagonal transport coefficients controlling thermoelectric 
effects (inset of Fig. 2). These protocols can be viewed 
as time-dependent analogues of the Seebeck and Peltier 
effects. We estimate the magnitude of the expected ef- 
fects in the simplest case of non-interacting fermions, for 
a ballistic or diffusive constriction with different geome- 
tries. We conclude that the proposed effects should be 
experimentally observable in the currently available set- 
up. 



The experimental realization of thermal and thermo- 
electric transport in cold atomic gases would open many 
possibilities. A particularly appealing opportunity of- 
fered by cold atomic gases is the study of these effects i) 
in the regime of high temperature where T is a sizeable 
fraction of the Fermi temperature Tp and ii) in the ab- 
sence of any phonon excitations. Purely electronic contri- 
butions to high-temperature transport and thermopower 
(as captured e.g. by Heike's formula) have been often 
discussed in the context of materials with strong elec- 
tronic correlations [21, 22] and are directly relevant to 
thermoelectric properties of oxide materials for example. 
In the solid state context however, separating the dif- 
ferent contributions of electrons and phonons is highly 
involved. A setup in which phonon contributions can 
be suppressed is therefore invaluable to reach a deeper 
fundamental understanding of thermoelectric transport 
of quantum interacting particles, which may well in turn 
provide guidance for better solid-state materials design. 
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Figure 1. Two-terminal transport setup. The left (L) and 
right (R) reservoirs are characterized by their temperatures 
Tl,r, chemical potentials /j.l,r and their thermodynamic coef- 
ficients k, Cfi, a, cf. Eq. (1). The constriction (region between 
dashed lines) is characterized by its matrix of linear response 
coefficients C. Due to the different temperatures and chemical 
potentials of the reservoirs, a particle (In) and entropy (Is) 
current flow in the constriction (in our conventions, I > 
corresponds to a flow from right to left). 

General framework. The transport setup under consid- 



eration is depicted in Fig. 1. Two reservoirs of fermionic 
atoms are connected by a constriction. The reservoirs 
are described by their temperature Tl,r and their chem- 
ical potential hl,r which determine the particle number 
N(T, n) and entropy S(T, /i) of each reservoir through the 
grand-canonical equation of state. Important thermody- 
namic properties of the reservoirs are the compressibility 
k, the dilatation coefficient a and the heat capacity C^.n 
(at constant fj, or N), defined by: 
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The particle and entropy currents flowing through the 
constriction I N = ±(N L -N R ), I s = £(S L -S R )[23] are 
related to the (small) chemical potential and temperature 
differences A/i = /i£ — [i R and AT = Tj, — Tr by: 



C is the matrix of transport coefficients associated with 
the constriction. Its symmetry is insured by Onsager's 
relations [24-26]. 

We consider a time-dependent process in which the 
reservoirs are prepared with given initial particle num- 
bers and temperatures, and equilibrate through exchange 
of particles and entropy through the constriction. 

Using (2) and the properties of the reservoirs given by 
(1), we derive equations ruling the time evolution of the 
particle and temperature imbalance: 
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These equations resemble the discharge of a capaci- 
tor (the reservoirs) in a resistor (the constriction), tak- 
ing into account thermal transport as well. The char- 
acteristic time-scale tq = k/Ch corresponds to the 
time-scale RC in a capacitor, with R = l/£n the re- 
sistance of the constriction and C (~ k) the capaci- 
tance [20]. The Lorenz number [27] of the constriction 



L = £ 22 /Al - (A2/A1) 2 is the ratio R/(TR T ) of the 
resistance to thermal resistance. The thermodynamic 
coefficient £ = C^/kT — (a/ft) 2 — Cn/kT character- 
izes the reservoir and has a form similar to L. Finally, 
S = ch/k — C12/C11 is the effective thermoelectric (See- 
beck) coefficient. Note that S has the dimension of the 
Boltzmann constant k R , while L and £ have dimension 
kg and k the dimension of the inverse of an energy. 

Two remarks are in order, (i) In the absence of any 
thermoelectric effect (S = 0), the time constants for 
particle and thermal relaxation are To and tqL/£, re- 
spectively. At low temperature, the usual form of the 
Wiedemann-Franz law, when applicable, dictates that 
L — > 7r 2 /3 and the thermodynamics of a free Fermi gas 
yields the same limit for £ — > 7r 2 /3. Hence, in this case, 
the Wiedemann-Franz law expresses the fact that the 
timescales for particle and heat relaxation are identical 
when S — 0, an interpretation which to our knowledge 
has not been formulated before. 

(ii) The coupling between heat and particle transport 
(off-diagonal elements of A) is determined by the effec- 
tive Seebeck coefficient S. In the present context, this 
coefficient has two competing contributions: one from 
the constriction S c = — Cy^l C\\ and one from the reser- 
voirs S r — a. jK. In particular, the presence of S r induces 
a coupling between thermal and electric transport even 
when the off-diagonal transport coefficient £12 of the con- 
striction can be neglected. 

On Fig 3, we display S C1 S r , L and £ as a function of 
T/Tp-, calculated within a Landauer-Biittiker [28-30] for- 
malism for a diffusive constriction (see below for details) . 
The plot illustrates the competition between the contri- 
bution of the reservoir S r and that of the constriction 
S c . In the case of Fig. 3, this results in a positive (neg- 
ative) total Seebeck coefficient in the absence (presence) 
of transverse harmonic confinement in the constriction. 
Note also on Fig. 3 the deviations from the Wiedemann- 
Franz law at finite temperature. 

The general solution of Eq. (3) reads, given an initial 
particle and temperature difference AJVq and ATq : 
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The inverse time-scales t ± = r A± are given by the 
eigenvalues of the transport matrix A 




In principle, eqs. (4) and (5) enable one to extract the 



thermodynamic and transport coefficients from experi- 
mental measurements. In particular, the Wiedemann- 
Franz law and its possible violations could be tested. 
However, thermoelectric effects are more difficult to ex- 
tract. This is due to the fact that in the presence of 
both a particle and temperature imbalance at t — 0, the 



time-evolution is typically dominated by the exponential 
decay involving the symmetric combination of exponen- 
tials, and terms responsible for offdiagonal transport are 
masked. 

Hence, we propose two specific protocols in order to 
reveal thermoelectric effects, estimate the expected mag- 
nitude of the signal and confirm experimental feasibility. 

Experimental proposal for off-diagonal (thermoelectric) 
transport. The first protocol focuses on the particle cur- 
rent induced by a temperature difference, and is a tran- 
sient analogue of the Seebeck effect. The system is pre- 
pared with equal number of particles in the two reser- 
voirs ANq — 0, but with a temperature difference (see 
inset of Fig. 2). The temperature difference can for ex- 
ample be prepared by heating one of the reservoirs (eg 
by using laser light) with a closed constriction, which 
prevents particle transfer. After reopening the constric- 
tion, the off-diagonal coupling between particle and heat 
transport will lead to an evolution of the particle number 
difference given by : 

AN(t) _ ST f k e-*/ T - - e-*/ T + AT 
N ~~ N A+ - A_ Tp 7 ' ( ' 

where No, Tp refer to reservoirs at equilibrium. Hence, 
a transient particle imbalance during equilibration is a 
'smoking gun' observation revealing the existence of a 
non-zero Seebeck coefficient. Since typically r + < t_, the 
particle imbalance first reaches an extremum at a time 
tmax — To(lnA+ — lnA_)/(A+ — A_) and then at long 
times relaxes exponentially to zero with a characteristic 
time r_. This behaviour is exemplified in Fig. 2. The 
sign of the particle imbalance is given by the sign of the 
effective Seebeck coefficient and thus depends on the con- 
sidered situation. If the reservoir effect is the dominant 
mechanism (S > 0), then particles tend at first to flow 
from the cold to the hot reservoir (see Fig. 2), due to its 
lower chemical potential. This behaviour is in contrast to 
the classical intuition of particles flowing from the high 
pressure (warmer) to the low pressure (colder) side. On 
the contrary, if the thermoelectric properties of the con- 
striction dominate (S < 0), then particles first flow from 
the hot to the cold reservoir. In both cases the tempera- 
ture imbalance equilibrates monotonically. As expected, 
the entropy current (Is = SIn — AT/RtT) is flowing 
from the hot to the cold reservoir, because the contribu- 
tion of thermal diffusion to entropy flow (second term) 
always dominates over the thermoelectric contribution 
from particle transport (first term). A second comple- 
mentary protocol can be proposed in which the system is 
prepared initially with a particle number imbalance but 
equal temperatures in the two reservoirs. This is the ana- 
logue of a Peltier experiment, in which temperature and 
particle number imbalances have reversed roles. 

In order to estimate the magnitude of the expected 
effects and assess whether they can be observed exper- 
imentally, we have considered the simplest case of non- 
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Figure 2. Inset: Sketch of the proposed protocol for detect- 
ing a transient thermoelectric Seebeck-like effect. The two 
reservoirs are prepared with equal particle numbers but with 
an initial temperature imbalance ATb. A transient particle 
imbalance signals a non-zero Seebeck coefficient. The arrows 
indicate the direction of the particle flow at short time in the 
different cases (see text). Mam: Particle and temperature im- 
balance (plain and dashed curves, respectively) as a function 
of time, for a 3D diffusive constriction with (blue) or without 
(red) transverse harmonic contribution (_L trap). The results 
depicted are for ATq/Tf = 25%, at an initial temperature 
T/T F = 0.4. 

interacting fermions and computed the transport coeffi- 
cients of the constriction using a Landauer-Biittiker for- 
malism, in the spirit of mesoscopic physics [28-30]. In 
this framework, the transport coefficients £y of the con- 
striction can be expressed in terms of the moments M n of 
an energy-dependent (dimensionless) transport function 
$(e) as: C n = 2M /h, C 22 = 2k%/h ■ M 2 /(k B T) 2 and 
£ 12 = 2k B /h ■ M 1 /(k B T), with: 

M n = jT°de*(e)(-|Q (e-/i)" (8) 

where / is the equilibrium Fermi function. The pre- 
cise form of the transport function $(e) depends on the 
dispersion relation and on the scattering mechanism in 
the constriction. We used a relaxation time approxima- 
tion with an energy-independent scattering time, and 
considered either a ballistic or a diffusive constriction 
(both cases can be realized experimentally [20] using 
speckle noise of adjustable strength). The constriction 
is three-dimensional, with or without harmonic confine- 
ment in the transverse direction, and the reservoirs are 
taken as three dimensionally harmonically trapped non- 
interacting Fermi gases. 

In order to observe the particle imbalance which is 
induced by the temperature imbalance, its amplitude 
should be sizable. For the typical parameters used in 
Fig. 2, The amplitude of the maximum AN(t max )/No 
is seen to exceed 10% , which makes it within reach of 
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Figure 3. Transport coefficients as a function of T/Tf for a 3D 
diffusive constriction, with (biue) or without (red) transverse 
confinement, and harmonic trapping in the reservoirs. 
(A)- Contributions from the reservoirs (SV, in units of ks) 
and from the constriction [S c , in units of fcs) to the total 
effective Seebeck coefficient S = S r + S c - (B)- Lorenz number 
L (in units of k%) of the constriction, and thermodynamic 
coefficient I (in units of k%) of the trapped reservoirs. 



current experimental setups [20]. To further quantify 
the effects, we define a parameter 77 = (AN (t max ) / N ) ■ 
(Tf/ATq) which can be understood as a thermoelectric 
efficiency of the system. It compares the maximum rela- 
tive particle imbalance to the initial temperature imbal- 
ance in units of the Fermi temperature. A similar quan- 
tity can be defined for the Peltier-like protocol, leading 
to an efficiency r]p = rj/l (hence expected to be smaller 
than 77 since I > 1 - see Fig. 3). Figure 4 displays 77 
as a function of T/Tp, comparing four cases - ballistic 
and diffusive and with or without transverse confine- 
ment in the constriction. In particular, one sees that 
77 ~ 50% at a commonly reached experimental temper- 
ature T/Tp — 0.4 [20], meaning that an initial relative 
temperature imbalance AT$/Tp = 20% with reservoirs 
containing iVo = 100000 particles would lead to a max- 
imal particle imbalance of the order of 10000 particles, 
which is experimentally sizeable. We note (Fig 4) that 
a very peculiar situation arises for a three-dimensional 
ballistic constriction with harmonic confinement. In this 
case, S r and S c perfectly compensate each other, leading 
to S = i] = 0. This is due to the fact that, in that 
case (assuming an energy-independent collision time), 
the transport function $(e) has exactly the same en- 
ergy dependence than that of the density of states in the 
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Figure 4. Efficiency 77 (see text) as a function of T/Tf for a 
ballistic or diffusive 3D constriction, with or without trans- 
verse confinement (with harmonic trapping of the reservoirs) . 



reservoirs (~ e 2 ). While this perfect cancellation relies 
on specific assumptions of our simplified modelling, the 
qualitative conclusion is expected to be robust: in this 
particular case the thermoelectric effects are expected to 
be quite small. 

To summarize, we have shown that thermoelectric ef- 
fects can be measured in cold atomic gases within the 
setup of Ref. [20]. In the current experimental tempera- 
ture regime, offdiagonal transport properties arise from a 
combination of the thermoelectric properties of the con- 
striction and the finite dilatation coefficient of the Fermi 
gas in the reservoirs. Fundamental questions about high- 
temperature transport can be addressed in this frame- 
work. 

Inter-particle interactions have been neglected here, 
implicitly assuming that they have been turned off us- 
ing a Feschbach resonance. A natural extension of this 
work is thus to include interactions, and to go beyond the 
simple Drude description of transport in the constriction. 
One can also think of reintroducing phonons in a con- 
trolled manner by simulating their action via a bosonic 
bath, an important ingredient for establishing contact 
with thermoelectric properties of solid-state materials. 
Other possible directions include the effect of a lattice 
in the constriction, or a modification of the geometry of 
the constriction, which would provide a fruitful relation 
to recent developments on thermoelectric devices in the 
mesoscopic physics context [31, 32]. 

Note added : As the writing of manuscript was being 
completed, we became aware of a recent preprint by H. 
Kim and D. Huse [33] considering spin and heat transport 
in a cold Fermi gas. 
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